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ABSTRACT 

The objects known as anomalous X-ray pulsars and soft gamma repeaters are 
commonly identified with magnetars, neutron stars with ultrastrong magnetic fields. 
The rotational history of these objects has, so far, revealed no evidence of free preces- 
sion. At the same time these objects do not generally appear to have magnetic axes 
nearly parallel or orthogonal to their spin axes. In this paper we show that the com- 
bination of these two observations, together with simple rigid-body dynamics, leads 
to non-trivial predictions about the interior properties of magnetars: either (i) elastic 
stresses in magnetar crusts are close to the theoretical upper limit above which the 
crustal matter yields or (ii) there is a "pinned" superfluid component in the magnetar 
interior. As a potentially observable consequence of these ideas we point out that, in 
the case of no pinned superfluidity, magnetars of stronger magnetic field strength than 
those currently observed would have to be nearly aligned/orthogonal rotators. 
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1 INTRODUCTION 

The magnetars are a subset of the neutron star popula- 
' tion with ultrastrong magnetic fields, typically greater than 
10 14 G, and long spin periods, typically 1 — 10 seconds. 
They manifest themselves observationally as anomalous X- 
ray pulsars (AXPs) and so ft gamma repeaters (SGRs); see 
IWoods fc Thompson! (120041 ) for a review. The latter class of 
objects are of particular interest, as periodicities have been 
observed in the tails of bursts from some objects, possibly 
providing the first ever evidence for excitation of neutron 
star normal modes. The theory of magnetars as very strongly 
magnetised neut ron stars was develop e d by Thompson & 
Duncan (see e.g. iThompson fc Duncanl (|l995. 1996. 200 il l), 
in which the decaying magnetic field powered the outbursts. 

In this paper we will concern ourselves with something 
that the magnetars don't do — they do not seem to undergo 
free precession. More precisely, we will base our analysis on 
the following two working assumptions, both motivated di- 
rectly by observations: (i) magnetars do not precess, and (ii) 
their magnetic pulsation axes are neither aligned nor orthog- 
onal to the spin axes. On the basis of these two observational 
results and simple rigid body dynamics we can then deduce 
that either (a) the stellar crust is very highly strained, right 
at or close to the limit predicted theoretically by detailed 
modelling, or (b) the magnetars contain a pinned superfluid 
component. We believe that either conclusion is interesting. 

The first assumption, the lack of precession, follows 



from timing studies. Precession would result in a periodic 
modulation in the time of ar rival (and possibly also sh ape) 
of the magnetar pulsations jjones fc A ndcrsson 200lj). In- 
deed, early observations of timing irregularities prompted 
Melatos to suggest precession was a generic feature of 
the AXP population (jMelatos 19 99|) . However, further ob- 
serva t ion of the AXPs ( Kaspi. Chakrabartv fc Steinbergerl 
Il999l: iKaspi et aljliobll : iGavriil fc Kaspil |2002| ) and SGRs 
( Woo ds et al . 2000) has ruled out precession at an amplitude 
level above the rms amplitude of magnetar timing noise. 
Small-amplitude precession, 'buried' inside the timing noise, 
is still a possibility and future timing data analysis should 
attempt to address this issue. (A strict lack of precession 
is not essential for our arguments: as will become clear, a 
small precession amplitude would imply that the magnetar 
spins about an axis close to a principal axis or the superfluid 
pinning axis; this is the essential point for our modelling). 

Our second assumption concerns the magnetar incli- 
nation angles, i.e. the angle 9b between their spin axes 
and magnetic axes; we assume that magnetars are neither 
aligned (9 b ~ 0) nor orthogonal (9b ~ ""/2). By talking 
of a 'magnetic axis', we are implicitly assuming there ex- 
ists a well defined axis about which the internal magnetic 
field is approximately symmetric; numerical simulat ions in- 
dicate this is likely to be the case (|Braithwaitell2009l ). Then, 
to make any observationally-motivated statement about the 
location of this axis, we also assume that this axis coincides 
with the axis defined by the magnetar pulsations, which is 
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probably a reasonable assumption given the known impor- 
tance of the ma gnetic field in determining th e surface tem- 
perature profile (jOzel. Psaltis fe K aspi 200l|). 

Given these assumptions, the result that the magnetars 
are not aligned rotators follows trivially from the fact that 
we see magnetar pulsations at all. That the two axes are 
not orthogonal can be argued on the basis of the harmonic 
content of the pulse profiles. An orthogonal rotator with two 
nearly identical anti-podal hot spots would, regardless of the 
observer's location, display a strong second harmonic in its 
pulse power spectrum, something which has not been re- 
ported t o be generic (but see th e data for IE 2259.1+586 in 
Fig. 2 of IGavriil fc Kaspil l|2002h . which does display strong 
harmonic structure). It therefore seems that the majority of 
the magnetars are neither nearly aligned nor nearly orthog- 
onal. However, one caveat needs to be added: the argument 
above, while safe for, say, a radio pulsar, where the emission 
comes from high magnetosp heric altitudes, is less secure for 
a magnetar. As described in lOzel. Psaltis fc Kaspil l|200lf ). if 
the magnetar radiation is produced at the stellar surface it- 
self, the combination of beaming and relativistic light bend- 
ing can result in even a single hot spot profile giving rise 
to a pulse with multiple peaks, making statements concern- 
ing the inclination angle very difficult to make. We regard 
the assumption of non-orthogonality between the spin and 
magnetic axes as the weakest link in our argument. 

Ours is not the first analysis of the precessio nal dynam- 
ics of magnetars. As noted above. iMelatod (|l999h attempted 
to interpret the irregular spin-down of AXPs as free preces- 
sion, but sustained observations did not confirm his model 
llKaspi. Chakrab arty fc Steinbergeill 19991 : iKaspi et al.ll200ll : 
IGavriil fc~K aspi 20o 3). M ore closely related to our work is 
that of lWassermanl (12003 ), who applied the same techniques 
as us (albeit without accounting for the presence of a super- 
fluid component) to work out conditions under which radio 
pulsars would not precess, and used these to made state- 
ments about the conditions under which precession is likely, 
and how the star's shape may change plastically over long 
timescales. As outlined above, the problem we address in 
this paper is different and perhaps simpler. The observa- 
tions indicate that the magnetars do not precess, which we 
then translate into constraints on the stellar structure. 



2 BASIC MODEL 

We will begin by assembling the most basic model that il- 
lustrates the key features of our argument. We assume a 
two-component neutron star model, consisting of a "crust" 
and a pinned neutron superfluid. In particular, the "crust" 
component includes the actual elastic crust, the charged 
fluids in the core as well as the portion of the neutron 
fluid not locked into the pinned component (vortex mutual 
friction would bring the non-pinne d superfluid into coro- 
nation with the charged component llAlpar. Langer fc Sauls! 
ll988l ; lAndersson. Siderv fc Comer! l20od )1. The two compo- 
nents are allowed to spin with different angular velocities 0,1 
and Ql. The former can be identified with the observed spin 
frequency fi* = fin^ (a "hat" denotes a unit vector). The 
crust's moment of inertia tensor can then be written 



This includes a spherical piece 7 C and the quadrupolar de- 
formations due to the elastic crust (A/ c ) and the interior 
magnetic field (AIb)- The two unit vectors hl,h B are as- 
sumed fixed in the crust frame, and we denote as /3 the 
angle between them. For the superfluid component we will 
assume a spherical moment of inertia tensor 



H 3 = U k 



(2) 



We will make the perfect pinning assumption, i.e. the pinned 
neutron vortices are physically immobilised with respect to 
the crust which in turn implies a superfluid angular velocity 
Ql — Q s nl which is fixed in the crust frame. 

It will be useful to parameterise the sizes of these con- 
tributions to the moment of inertia tensor. If 1$ — I c + I s 
denotes the spherical moment of inertia of the whole star, 
we can define 

AJ C AI B Is 

e c = —j—, 6b = — — , e s = -p. 
Jo Jo Jo 

The deformation e c is sou rced by strains in the solid 
crust. The recent calculation bv lHaskell. Jones fc Anderssonl 
(2006) place upper limits on this of 
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where we have parameterised the crustal breaking strain Ubr 
in term s of the value found in th e state-of-the-art calcula- 
tions of iHorowitz fc Kadaul (|200gh . 

For the magnetic deformation a somewhat crude es- 
tima te (but consistent w it h mo r e rigorous calcu l ations , 
e.g. lAkgun fc Wassermanl <|2008h : IColaiuda et"ail l|2008h : 
lLander fc Jones! (|2009l VHs 
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where B is a volume average of the internal magnetic field 
and H rss 10 15 G if the core sustains type II proton super- 
conductivity, otherwise H = B. The deformation ea can be 
either positive or negative depending on the relative strength 
between the poloidal and toroidal ma gnetic field components 
l|Mestel fc Takharlll972l : ICutlerll2002r i. 

Clearly, given the previous estimates, the magnetic and 
elastic deformations are small numbers; the same is not nec- 
essarily true for e s . If pinning occurs only in the inner crust, 
only a few percent of the stellar superfluid contributes to I s , 
giving e s « 10 -2 . However, if vortex pinning is efficient in 
the neutron star core due to the interaction with the mag- 
netic fluxtubes then e s could even be of order unity. Either 
way, e s is likely to be much greater than either of e c and eb- 

We will now easily find the non-precessional solutions. 
The total angular momentum of our star is given by the sum 
of the crustal and superfluid contributions 

j i = ii'n i + i?a , i . (6) 

Let 9 C and 8b denote the angles that the crustal and mag- 
netic deformations make with the angular velocity vector, 
so that n^nf = cos0 c and n B nf = cos 0b- The angular mo- 
mentum can then be written 

T = Q. (^IcTiQ + AI C cos 8 c nl + AIb cos^b^b^ + 



I C S V + AI c n c n 3 c + A/ B w B ng. 



(1) 



The Euler equations of rigid body dynamics then give 
j i + e ijk n j J k =0, 



(7) 



(8) 
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where the time derivative is evaluated in the rotating crust 
frame. For a non-precessional and therefore stationary solu- 
tion we must have J ! = which implies that J ! must be par- 
allel to i.e. there exists a scalar A such that J % = XQuq. 
We therefore obtain the relationship 

(I c - \)n}i + AI c cos8 c n c + AI B cos6 B n B + Is^-hl = (9) 

between the four unit vectors n c , Ab, fin, n a . 

Equation Q is crucial-we can easily derive our two key 
results from it. Firstly, consider the case where there is no 
pinned superfluid. Setting 7 S = we obtain 

^— — n n + e c cos 9 c fi c + e B cos 9 B n' B = 0. (10) 

This shows that the three unit vectors n c , n B , nn, are copla- 
nar. We can therefore make the decomposition into a pair 
of basis vectors h\ , h\ such that 



(11) 
(12) 
(13) 

Inserting into equation (|10[) and projecting along h\ leads 
to the non- precession condition 



n B — sm o B n 1 + cos p B n 3 , 
n c — sin 6 c h\ + cos 9 c h\, 

~ i » i. 



' sin 28 c = — e B sin 29 B . 



(14) 



Thus, in the case of zero pinning, non-precession implies a 
simple relation between the crustal an d magnetic deforma - 
tions. The same result was obtained by I Wassermanl |2003). 

Now consider the case where there is pinning. We de- 
fine the angle n'np = cos# s . Taking the vector product of 
equation ([9| with n n and squaring the result leads to 

\ 2 

2e s -^sin(9 s j = (ecsin20 c ) 2 + (e B sin20 B ) 2 

+ 8e c £B cos 6 C cos # B (cos /3 — cos 8 C cos 6 B ). (15) 

This is a rather cumbersome relation between the angles 
that the superfluid, the crustal deformation, and the mag- 
netic axis make with the rotation axis. However, from this we 
see one important thing: the star rotates about an axis very 
close to the superfluid pinning axis, with the angle of mis- 
alignment being of order of the larger of AI C /I B and AI B /I B . 
In the notation of Q, we therefore have 



(16) 



Note that in the pinning scenario, the magnetar non- 
precession does not place any constraint on the crustal de- 
formation. Also, it is possible for the star to spin about the 
crustal principal axis (as defined by eqn. (|T4j)) and have 
a pinned superfluid, provided the pinning axis lies exactly 
along the spin axis (9 S = 0). This would correspond to both 
our scenarios being realised. However, there is no apriori 
reason to expect the symmetry axis of the pinning sites to 
coincide with the spin axis, so this is a clearly special case. 

It is also interesting to note the well known result that 
the non-precessing state is the one that mimimis e s the ki- 
netic energy at fixed J (|Landau fc Lifshitzl 1 19761 ; IShaham] 
1 19771 ). Therefore, an isolated star, dissipating kinetic energy, 
would evolve to such a state on the dissipation timescale. 

Finally, note that this last point has repercus- 
sions for the so-called Mestel- Jones spin-flip mechanism 



jMestel fc Takharlll972i ; rJones!ll975l ; rCutleij|2002i ). in which 
a star with a dominantly prolate magnetic deformation of 
its inertia tensor (i.e. e B < 0) is driven to a configuration 
with the magnetic axis orthogonal to the spin axis. Clearly, 
the presence of a pinned superfluid completely suppresses 
this, assuming that the magnetic and pinning axes remain 
fixed with respect to each other (as in our model). Accord- 
ing to the above discussion, the minimum energy state for a 
star containing a pinned superfluid has 8 S ~ 0, i.e. the pin- 
ning axis must remain close to the spin axis, preventing the 
migration of the magnetic axis to an orthogonal location. 



3 EXTENSIONS TO THE BASIC MODEL 

The previous model can be extended in terms of realism by 
incorporating additional physical effects. These include (i) 
rotational "bulges" in the moments of inertia, and (ii) the 
action of a non-dissipative torque due to the external dip ole 
field (the so-called "anomalous" torque l| Goldreichl 1 1 9 70l ) V 

With the rotational deformations included the previous 
moment of inertia tensors become 

I? = h8 l] + AInnhnh + A7 c n c n c + AI B n B h 3 B , (17) 



I s 8 13 + AI S . 



(18) 



For a slowly rotating object like a magnetar we would expect 
AIn 7 C and A7 S I s . The total angular momentum is 

J 1 = n | (7 C + A7 n ) fiii + AI C cos (9 c n c 

+A7 B cos6> B n B } + (7 s + A7 s )a i n*. (19) 

This expression shows that the centrifugal deformations are 
not expected to play any significant role since they are sim- 
ply "absorbed" in the pre-existing spherical pieces. 

The second modification to the basic model, the anoma- 
lous torque A^n, will appear as a source term in the Euler 
equation ((H) . If the exterior space is assumed perfect vacuum 
and the magnetic field dipolar this torque can be written as 



A/a n = IofleA(nnmi)e ljk D.j7hk, 



(20) 



where m l is the external magnetic dipole moment and e\ is 
an effective "deformation" 



and R is the stellar radius. Expressed in terms of the field 
strength Bd on the magnetic pole we obtain the estimate 



IO" 7 



B d 



10 15 G 



(22) 



Clearly, the impact of the anomalous torque on the stel- 
lar dynamics is important only for magnetar-strength fields. 
This was the central ide a in the s o-calle d radiative preces- 
sion model advocated bv lMelatosI j 19991 ). 

With the inclusion of the anomalous torque the equa- 
tion of motion ((HJ takes the form 

j l + e l3k Q,j \j k - e A I tt(h n mi)rh k ^ = 0. (23) 

Precession will not occur provided 

J 1 = XttfiQ + eAloQ(h Q mj)rn' . (24) 
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It is natural to assume that the interior and exterior mag- 
netic fields share the same symmetry axis (see discussion in 
the next section). Setting ra l = h B in (|24|l we recover the 
previous results (114[) - (|16| l with the recalibrated deformations 
£b — > e B - £A and e s — > e s + A7 s /7 . 



4 DISCUSSION 

First consider equation <|14[) . which applies in the limit of no 
pinning. The magnetic ellipticity cb is at least as large as 
would be estimated by substituting the external magnetic 
field strength into equation ((5j), i.e. for a typical magnetar 
e B ~ 1(T 6 for B = 10 15 G. In reality this number could 
be significantly larger if we take at face value the most re- 
cent calculations of stable magnetic equilibria in neutron 
star models which suggest that the inter ior field B could b e 
much stronger than the exterior field Bd <|Braithwaitell2009l ). 

That the magnetars do not seem to be either aligned or 
orthogonal rotators means that the angle 9 b is neither close 
to zero nor close to n/2, and so the trigonometric factor on 
the right hand side of (|14|l is of order unity. It then follows 
that the crustal deformation must be of the order of the 
magnetic one: e c ~ e B > 10 -6 , or possibly even larger, if 9 C 
is small. The value e c ~ 10~ 6 is rather significant- it is close 
to the maximum allowed deformation predicted by detailed 
modelling of the crust, i.e. equation (Q. 

It follows that, in the absence of superfluid pinning, 
magnetar crusts are maximally, or close to maximally, 
strained. This ide a sits well with the standard model o f 
magnetar activity l|Thompson fc Duncan|[l995. 1996. 200lh 
where the energy budget of the flaring activity is supplied by 
the evolving magnetic field. The crust itself acts as a gate, 
with crustal cracking acting as a trigger for the magnetic 
reconnection events. The maximally strained crust that ap- 
plies in this scenario makes such events perfectly natural. 

One might worry that one could evade our conclusion 
of a strongly strained crust by relaxing the assumption that 
the internal and external magnetic fields are aligned (i.e. 
rn 1 7^ n B ). Indeed, if the internal magnetic axis lies close 
to the rotation axis (making the angle 9 B small), equations 
(|14p allow for a e c well below the maximum. However, a mis- 
alignment between the internal and external magnetic fields 
will generate tangential magnetic stresses, which can only be 
balanced by elastic ones. If both the internal and external 
fields are of strength ~ B, and if the field transitions from 
one geometry to the other over a shell of thickness AR < R, 
the tangential magnetic stresses will be the order of B 2 /AR. 
This stress is greater than the B 2 /R that the crust needs 
to balance in our original scenario with aligned internal and 
external fields. Clearly, relaxing this assumption does not 
allow one to avoid our conclusion of a highly strained crust. 
Additional evidence supporting the alignment between the 
interior and exterior magnetic fields comes from simulations 
of stable magnetic equilibria l|Braithwaitell2009l ). 

Now consider the case with superfluid pinning (equa- 
tions (|15[) and (|16|) ). In this case no constraint is placed 
on the crustal deformation. In fact, e c could be set to zero, 
with no qualitative change in the non-precessional solution. 
The only constraint on stellar parameters that the pin- 
ning scenario would imply is that the pinning force itself 
is strong enough to sustain the small misalignment between 



the crustal and superfluid ang ular velocity vectors given by 
equation (fl6|) . As discussed bv lLink fc Cutlerl (|2002l ) the rel- 
ative velocity between the superfluid and the crust gener- 
ated by this misalignment produces a Magnus force on the 
pinned vortices. Making a simple estimate, the relative ve- 
locity generated by the misalignment is of order An ~ QR9 S . 
Parameterising with magnetars in mind gives 

This is to be compared with the critical velocity f or un- 
pinnin g. The estimates of this have been collected in I Jonesl 
(2010) who finds a lower unpinning threshold ~ 10 4 cm/s S> 
Av. Thus, the finite strength of vortex pinning is no obstacle 
to building our non-precessing magnetar models. Further- 
more, the weak Av also ensures that the non- precession con- 
fi guration is immune to the superflu i d inst ability discussed 
in lGlampedakis. Andersson fc Jonesl (|2009l ). 

There is one special case that deserves some comment: 
when the crustal and magnetic axes are aligned, so that 
n l c = n B . In this case the star is biaxial rather than triaxial, 
and precession would correspond to the magnetic axis rotat- 
ing in a cone of half-angle « 9b about the fixed vector J 1 , 
with a rotation about h B at the precession frequency super- 
imposed. As is well known, such a precession would h ave no 
signature in the pulsations (I Jones fc Anderss on 2001). Such 
a seemingly non-precessing star could have no pinned super- 
fluid and any amount of crustal strain (even zero-in this case 
the star is still biaxial), and so would not fit into the scheme 
described above. How can we be sure that magnetars do not 
fall into this special case? Firstly, the absence of precessional 
modulation in the pulsations would require the hot-spot it- 
self to be exactly axisymmetric about the magnetic axis. 
Secondly, it would be surprising if the crustal strains were 
symmetric about the magnetic axes, as the crustal strain 
field is likely to have retained some memory of a more oblate 
shape inherited from earlier in its life, before it was spun- 
down to the long periods of the magnetars. Thirdly, such a 
large-amplitude precession woul d be expected to damp du e 
to internal dissipative processes jjones fc Andersson| [2001). 
We therefore feel it is very unlikely that this special case is 
common in the magnetar population. 

Having arrived at two different non-precession possibil- 
ities, described by equations (|14fl and (|15[) , it is natural to 
ask the question of how one might hope to distinguish be- 
tween them. As already noted, the no pinning case, with 
its maximally strained crust, sits naturally with the mag- 
netar flare model, as quake events in the crust are required 
to trigger magnetic reconnection events. However, this is 
hardly decisive-diffusion of the magnetic field may, in either 
scenario, play a key role in triggering the bursts. 

Nevertheless, one intriguing possibility does suggest it- 
self. In the scenario with no superfluid pinning the crustal 
strain predicted just happens to lie at the upper end of what 
is believed to be possible on the basis of detailed crust mod- 
elling. Suppose this is no coincidence. It may be the case that 
there exist even more strongly magnetised stars; their mag- 
netic deformations would necessarily be larger than their 
crustal ones. By eqn. (fTl)) it follows that such stars, provided 
they are oblate (es > 0), would be nearly aligned rotators. 
There would be a selection effect against seeing such mag- 
netars. So, observation of a tail of strongly magnetised but 
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nearly aligned rotators would argue in favour of there being 
no superfluid pinning in magnetars. Similarly, if the mag- 
netic deformation were prolate (cb < 0), strongly magne- 
tised stars would be orthogonal rotators, although as noted 
above, the current observations do not suggest this to be the 
case. Note, however, regardless of the sign of ea, the absence 
of such aligned/orthogonal populations would not necessar- 
ily imply that superfluid pinning occurs in magnetars; Na- 
ture may simply not provide neutron stars with magnetic 
field strengths higher than those currently observed. 

The presence of a pinned superfluid is also related to 
the timing profile of magnetars. Indeed, several AXPs are 
known to have undergone one or more glitches with a frac- 
tional f2-jump compara ble to the Vela-type glitches in radio 
pulsars (for a review see lRim. Kaspi fe Gavriil |2008l )). The 
occurrence of large glitches by itself suggests the presence 
of a pinned superfluid reservoir somewhere in the stellar 
interior, capable of transfering angular momentum to the 
crust if some mechanism could cause vortex unpinning. So 
far no glitches have been detected in SGRs. This should 
not come as a surprise given that these objects, when in 
quiescence, are actually fainter than the AXPs as well as 
"louder" in terms of timing noise. As a result, the phase- 
coherent timing of SGRs required for the identification of 
glitches becomes a highly problematic procedure. We can, 
nevertheless, speculate on physical grounds as to why SGRs 
may not glitch. One possibility could be the absence of vor- 
tex pinning in the core, if the core is the only region where 
pinning could in principle occur. The interior magnetic field 
in SGRs could be sufficiently strong as to exceed the critical 
field H C 2 w 10 16 G above which proton su perconductivity is 
suppressed l|Bavm. Pethick fc Pineslll969h . The destruction 
of superconductivity would obviously imply the absence of 
magnetic fluxtubes on which the vortices could pin. If this 
scenario were to be true then the known SGRs would obey 
the non-precession condition ()14[) and would have maximally 
strained crusts. 

It is interesting to consider the impact that a hypo- 
thetical discovery of magnetar precession would have on 
our principal conclusions. As we discussed in the introduc- 
tion, current observations do not exclude the possibility that 
magnetars could be undergoing small-amplitude precession, 
masked by their timing noise. This precession could be either 
fast or slow (i.e. with a period, respectively, comparable to 
or much longer than the spin period) depending on whether 
there is a pinned superfluid. A detection of the former type 
of precession would amount to clear evidence of the existence 
of a pinned superfluid component within the star. Alterna- 
tively, detection of slow precession would mean that (i) there 
is no pinned superfluid, and (ii) that the star's rigid body 
dynamics deviates only slightly from the non-precession con- 
dition (|14p . if the precession amplitude is sufficiently small. 
Thus, even in this case we would be able to make the same 
predictions for the magnetar structure as before. 

To sum up, we have argued that the high magnetic 
strains that must necessarily exist in magnetars have inter- 
esting implications for their structure, telling us something 
about the level of crustal strain or the superfluid state of 
matter in their interior. Our modelling assumes a lack of 
significant free precession in the magnetars, and that their 
magnetic axes, as traced by the pulsation axes, are neither 
nearly aligned nor nearly orthogonal to their spin axes. We 



therefore conclude by pointing out the importance of better 
observational data to either confirm or reject these assump- 
tions. In particular, our work shows that there would be a 
great deal of interest in attempting to construct detailed 
models of the geometry of magnetar emission, despite the 
technical difficulties involved. There would also be interest in 
searching for an ultra-highly magnetised but nearly aligned 
or orthogonal tail to the magnetar population. Clearly, more 
detailed observations of magnetar timing and pulse geome- 
try could give us a unique insight into neutron star interiors. 
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